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ABSTRACT 

The Shapley Supercluster Survey is a multi-wavelength survey covering an area of ~23 
deg^ (~ 260 Mpc^ at z=0.048) around the supercluster core, including nine Abell and two 
poor clusters, having redshifts in the range 0.045-0.050. The survey aims to investigate the 
role of the cluster-scale mass assembly on the evolution of galaxies, mapping the effects of 
the environment from the cores of the clusters to their outskirts and along the filaments. The 
optical (ugri) imaging acquired with OmegaCAM on the VLT Survey Telescope is essential to 
achieve the project goals providing accurate multi-band photometry for the galaxy population 
down to m*-i-6. We describe the methodology adopted to construct the optical catalogues and 
to separate extended and point-like sources. The catalogues reach average 5cr limiting mag¬ 
nitudes within a 3" diameter aperture of Mgn=[24.4,24.6,24.1,23.3] and are 93% complete 
down to Mgn=[23.8,23.8,23.5,22.0]mag, corresponding to ~ m*r-i-8.5. The data are highly 
uniform in terms of observing conditions and all acquired with seeing less than 1.1 arcsec 
full width at half-maximum. The median seeing in r-band is 0.6 arcsec, corresponding to 
0.56kpch^y at z=0.048. While the observations in the u, g and r bands are still ongoing, the 
f-band observations have been completed, and we present the /-band catalogue over the whole 
survey area. The latter is released and it will be regularly updated, through the use of the Vir¬ 
tual Observatory tools. This includes 734,319 sources down to /=22.0mag and it is the first 
optical homogeneous catalogue at such a depth, covering the central region of the Shapley 
supercluster. 

Key words: methods: data analysis - methods: observational - catalogues - virtual obser¬ 
vatory tools - galaxies: clusters: general - galaxies: clusters: individual: A 3552 - galaxies: 
clusters: individual: A 3554 - galaxies: clusters: individual: A 3556 - galaxies: clusters: indi¬ 
vidual: A 3558 - galaxies: clusters: individual: A 3559 - galaxies: clusters: individual: A 3560 
- galaxies: clusters: individual: A 3562 - galaxies: clusters: individual: AS 0724 - galaxies: 
clusters: individual: AS 0726 - galaxies: clusters: individual: SC 1327-312 - galaxies: clus¬ 
ters: individual: SC 1327-313 - galaxies: photometry. 


1 INTRODUCTION 


E-mail: mercurio@na.astro.it 


The main aim of the Shapley Supercluster Survey (ShaSS) 
is to quantify the influence of hierarchical mass assembly on 
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galaxy evolution and to follow such evolution from filaments 
to cluster cores, identifying the primary location and mecha¬ 
nisms for the transformation of spirals into SOs and dEs. The 
most massive structures in the local Universe are superclusters, 
which are still collapsing with galaxy clusters and groups fre¬ 
quently interacting and merging, and where a significant num¬ 
ber of galaxies are encountering dense environments for the 
first time. The Shapley supercluster (hereafter SSC) was cho¬ 
sen becau se of i) the pecu l iar cl u ster, galaxy and bary o n over¬ 
densi ti es I Scaramella et all 19^; Ra^chaufiiu^ 19891: Fabiar 


1991 


iRavchaudhurv et ~ Il99ll : lOe Filippis. Schindler & Erber 


200511 : ii) the relative dynamical immaturity of this super¬ 


cluster and the possible presence of infalling dark matter 
haloes as well as evidence of cluster-cluster mergers (e .g. 
BMdelli_et_^ I|1994 lOuintana et al.lll995l: iBardel l i et~^ll998alltt 


Kulj_&_Bdhrin 2 et 1^99: guintanaj^am^o & Rei ^neggeil 200C : 
Bardelli. Zucca & Baldi200ll ' lDrinkwater et al.l2004 ): iii) the pos¬ 


sibility that it is the most massive bound structur e known in the Uni¬ 
verse , at least in the lOMpc central region (see lPearson & Batuskil 
l2013h . These characteristics make the SSC an ideal laboratory for 
studying the impact of hierarchical cluster assembly on galaxy 
evolution and to sample different environments (groups, filaments, 
clusters). Furthermore, the redshift range of this structur e (0.033< 
z <0.060 lOuintana'^ ^1 19951 1 19971 : IProust et^ 1 20061) makes it 
feasible to measure the properties of member galaxies down to the 
dwarf regime, providing that the observations reach the suitable 
depth. A detailed discussi on of the scientific aspects of the survey 
is given in Merluzzi et al. I 20 TI) . 

Although, the SSC has been inves tigated by numerous 
authors since its disco very (Shapley 193(lh bot h for its cosmo¬ 
logica l implications 1 Scarameila_et_alJ [l989|; Ravchaudhury l 

!9a; 


19891: [Plionis & Valdamin? 


1991 


lOuintana et al.l 1 1995 


Kocevski, Mullis & Ebelingl I 2 OO 4 ; Feindt et 20131 and ref- 


Kull & Bohringeil Il999l: 

Bardelli et alJ I 2 OOOI: iFinoguenov et al.l 

2004; Rossetti et al.l 2005 

; [Munoz & LoebI 2008). none of them 


could systematically tackle the issue of galaxy evolution in the 
supercluster environment due to the lack of accurate and homo¬ 
geneous multi-band imaging covering such an extended structure. 
ShaSS aims to fill this gap measuring the integrated [magnitudes, 
colours, star formation rates (SFRs)] and internal (morphological 
features, internal colour gradients) properties of the supercluster 
galaxies. 

ShaSS will map a region of ~ 260 h^^ Mpc^ (at z=0.048), cen¬ 
tred on the SSC core, which is constit uted by three Abell clusters: 
A 3558 (z=0.048, Melnick & Ouintana ll98l[ Metcalfe, Godwin & 
SDenser lT987l : Abell richness R=A, Abell, Corwin & 01owin lT98^ . 
A 3562 (z=0.049, R=2) and A 3556 (z=0.0479, R=0); and two poor 
clusters SC 1327-312 and SC 1329-313. The present survey cov¬ 
ers also six other Abell clusters: A 3552, A 3554, A 3559, A 3560, 
AS 0724, AS 0726, as shown in Fig.[T] The survey boundaries are 
chosen not only to cover all 11 clusters and the likely connecting 
filaments, but also to extend into the fiel d and to map the struct ures 
directly connected to the SSC core. In iMerluzzi et al.l ( l2015h we 
derived the stellar mass density distribution based on supercluster 
members showing that all the clusters in the ShaSS area are em¬ 
bedded in a common network and identified a filament connecting 
the SSC core and the cluster A 3559 as well as the less pronounced 
overdensity extending from the SSC core towards A 3560. 

The data set of the survey includes optical {ugri) and NIR 
{K) imaging acquired with (VST) and Visible and Infrared Sur¬ 
vey Telescope for Astronomy (VISTA) respectively, and opti¬ 


cal spectroscopy with AAOmega. At present the i-band imag¬ 
ing and AAOmega spectroscopic surveys are completed, while 
the other observations are ongoing. In addition, the recent pub¬ 
lic relea se of data from the Wide-field Infrared Survey Explorer 
(WISE, IWright et al.l I 2 OI 0 I) provides photometry at both near- 
IR (3.4,4.6/rm) and mid-IR (12,22yum) wavelengths, allowing in¬ 
dependent measurements of stellar masses down to AI=10®Mo 
at lOcr and SFR down to 0.46Moyr“* at lOcr (0.2Moyr“* at 
5cr). Finally, in the 2-3 deg^ of the SSC core, panoramic imag¬ 
ing in the UV (Galaxy Evolution Explorer, GALEX), optical 
(ESO Wide Field Imager, WFI), NIR (UKIR T/WFCAM) and mid- 
infrared (Spitzer/MIPS) are also availa ble jMercurio et alJlT006l : 
IMerluzzi et alj2oT^ : lHaines et al.ll201lh . 

The optical survey, whose coverage is indicated by the 1 deg^ 
boxes in Fig. [T] will enable us to i) derive accurate morphologies, 
structural parameters (Slogr^ ~ 0.04 and Sns„- ~ 1) as well as 
detect some of the observational signatures related to the differ¬ 
ent processes experienced by supercluster galaxies (e.g. extrapla- 
nar material); ii) estimate acc urate colours, photo-zs (dz < 0.03, 
see lChristodoulou et al.ll2oT2l) and stellar masses; iii) evaluate the 
SFRs and resolve the star forming regions at least for the subsample 
of brighter galaxies. 


The survey depth enables global and internal physical proper¬ 
ties of Shapley galaxies to be derived down to m*+6. In the first 
case of obtaining accurate measurements of aperture photometry 
and colours, we require signal-to-noise ratios (SNR) of 20 in all 
four bands for SSC galaxies down to m*+6. Secondly, for the mor¬ 
phological analysis and resolving internal properties and structures 
there is a mo re stringent requirement of SNR~ 100 (in a 3" diameter 
apertu re, see IConselice. Bershadv & Jangren|[2000l : lHiussler et ^ 
l2007li for the deeper r-band imaging. For this reason we are col¬ 
lecting the r-band imaging under the best observing conditions, 
with a full width at half-maximum (FWHM)~0.8arcsec or better, 
corresponding to 0.75 kpc at z=0.048. Additionally, the r imaging 
is fundamental to our weak lensing analysis, to ensure a sufficient 
density of lensed background galaxies with shape measurements. 


With these data it will be possible to separate the different 
morphological types, trace ongoing SF, reveal recent interaction or 
merging activities and thus obtain a census of galaxies whose struc¬ 
ture appears disturbed by the environment (e.g. Sc^lataet^ .I 2 OO 7 I: 


Munoz;Mateoset_^ , 

20l4lLotz et alll2008L 


2004 iKleiner et alJ I20l4 Holwerda et al.l 
201 iL and references therein). 


In order to achieve the scientific goals of the survey, accurate 
photometry is required. This implies a clean source catalogue con¬ 
taining, together with the measured photometric properties, indica¬ 
tors of the reliability of these measurements. This paper describes 
the methodology used to produce the photometric catalogues and 
the adopted procedures. 

Observations are overviewed in Sect. In Sect. [3 we de¬ 
scribe the construction of the catalogues, the criteria to classify 
spurious objects and unreliable detections, the flags adopted in the 
catalogues and the procedure for star/galaxy separation. The ac¬ 
curacy and completeness of the derived photometry is discussed in 
Sect.|4] Each parameter of the released i-band catalogues is detailed 
in Sect.|5]and the summary is given in Sect.|6] 


Throughout the paper, we assume a cosmology with f2n,= 0.3, 
flA= 0.7 and Ho= 70 km s“’ Mpc’. According to this cosmological 
model 1 arcmin corresponds to 56.46 kpc at z=0.048. The magni¬ 
tudes are given in the AB photometric system. 
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Figure 1. VST fields mapping the ShaSS region. Dots indicate the super¬ 
cluster members in the range 13500-16000 km s*, taken from 

literature. Note that the literature redshift coverage is not uniform across the 
region shown. Black crosses show the cluster centres. 

2 VST OBSERVATIONS AND DATA REDUCTION 

The optical survey (PI: P. Merluzzi) is being carried out using 
the Italian INAF Guaranteed Time of Observations (GTO) with 
OmegaCAM on the 2.6m ESO VLT Survey Telescope located at 
Cerro Paranal (Chile). The Camera has a corrected field of view 
of l°xl°, corresponding to ~ 3.4x3.4 h^g Mpc^ at the supercluster 
redshift, sampled at 0.21 arcsec per pixel with a 16kxl6k detector 
mosaic of 32 CCDs, with gaps of 25-85 arcsec in between chipsQ 

Each field is observed in four bands; ugri. To achieve the re¬ 
quired depth, total exposure times for each pointing are 2955s in u, 
1400s in g, 2664s in r and 1000s in i. To bridge the gaps a diagonal 
dither pattern of five exposures in u, g and i, and nine exposures 
in r is performed, with step size of 25" (15" for the r band) in 
X, and 85" (45" in r) in the Y-direction. The total area is covered 
by 23 contiguous VST pointings overlapping by ~3' as shown in 
Fig.m where dots denote the spectroscopic supercluster members 
(13500 < Vi, < 16000 km s“*) available from literature at the time 
of the survey planning. The X-ray centres are indicated by crosses 
for all the known clusters except AS 0726, whose centre is derived 
by a dynamical analysis. 

The survey started in 2012 Eebruary and will be completed in 
2015 (spanning ESO periods P88-P95), provided that all the fore¬ 
seen observations are carried out. At present, the survey coverage 
differs for each band with only the i-band observations available for 
the whole area. This implies that the results concerning the quality 
of the photometry (depth, completeness, accuracy) are based on a 
representative subsample of the final catalogues for ugr (48%, 43% 
and 61%, respectively) bands and for the whole catalogue in i band. 

The data are collected on clear and photometric nights with 
good and uniform seeing conditions. In Fig. [2] we plot the see¬ 



Seeing (arcsec) Seeing (arcsec) 


Figure 2. Seeing distribution in the four observed bands (see text). 


ing values of 11, 10, 14 and 23 fields in ugri, respectively. Out 
of the observed fields about 80% (gi) and 90% (r) are acquired 
with FWHM< 0.8 arcsec, with the median seeing in r band equal 
to 0.6 arcsec, corresponding to 0.56 kpc h^g at z=0.048. The u band 
is characterized by a slightly poorer seeing. We discuss the effect 
of the seeing on the aperture magnitudes in Sect. 13.1.^ 

_ ^The data reduction is already described in Merluzzi et al. 

I 2 OI 5 I Here we summarize the main steps for the reader’s conve¬ 
nience. 

Images are reduced a nd com bined using the VST-Tube imag¬ 
ing pipeline (Grado et al. l2012h . developed ad hoc for the VST 
data. The pipeline follows the standard procedures for bias sub¬ 
traction and flat-field correction. A normalized combination of the 
dome and twilight flats, in which the twilight flat is passed through 
a low-pass filter first, were used to create the master flat. A gain har¬ 
monization procedure has been applied, finding the relative CCD 
gain coefficients which minimizes the background level differences 
in adjacent CCDs. A further correction is applied to account for the 
light scattered by the telescope and instrumental baffling. This is 
an additional component to the background, which, if not corrected 
for, causes a position-dependent bias in the photometric measure¬ 
ments. This component is subtracted through the determination and 
the application of the illumination correction (IC) map. The IC map 
is determined by comparing the magnitudes of photometric stan¬ 
dard fields with the corresponding SDSS-DR8 (Sloan Digital Sky 
Survey-Data Release 8) PSF magnitudes. 

For the i band a correction is required because of the fringe 
pattern due to thin-film interference effects in the detector from sky 
emission lines. The fringing pattern is estimated as the ratio be¬ 
tween the Super-Flat and the twilight sky flat, where Super-Flat is 
obtained by overscan and bias correcting a sigma clipped combi¬ 
nation of science images. The fringe pattern is subtracted from the 
image, applying a scale factor which minimizes the absolute differ¬ 
ence between the peak and valley values (maximum and minimum 
in the image background) in the fringe corrected image. 


* More details on the camera are available at 
http://www.eso.org/sci/facilities/paranal/instruments/omegacam/instHlDentPhotometric calibration onto the SDSS photometric sys¬ 
tem is performed in two steps: first a relative photometric cali- 
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Band 

Colour term 

Colour 

Extinction 

Zero Point 




Coefficient 


M 

r 

c 

A 

ZP 

u 

0.026±0.019 

u-g 

0.538 

23.261±0.028 

g 

0.024±0.006 

g-i 

0.180 

24.843±0.006 

r 

0.045±0.019 

r-i 

0.100 

24.608±0.007 

i 

0.003±0.008 

g-i 

0.043 

24.089±0.010 


Table 1. Absolute photometric calibration coefficients. 


bration among the exposures contributing to the final mosaic im¬ 
age is obtained through the comparison of the magnitudes of the 
bright unsaturated sta rs in t he different exposures, using the soft¬ 
ware SCAMP (Bertin 1200^^1 : then the absolute photometric cali¬ 
bration is computed on the photometric nights comparing the ob¬ 
served magnitude of stars in photometric standard fields with SDSS 
photometry. For those fields observed on clear nights, we take ad¬ 
vantage of the sample of bright unsaturated stars in the overlap¬ 
ping region between clear and photometric pointings and by using 
SCAMP, each exposure of the clear fields is calibrated on to the 
contiguous photometrically calibrated field. The magnitude is then 
calibrated by adopting the following relation: 

M'= M + yC + AX + ZP, (1) 

where M is the magnitude of the star in the standard system, M' is 
the instrumental magnitude, y is the coefficient of the colour term, 
C is the colour of the star in the standard system, A is the extinction 
coefficient, X is the airmass and ZP is the zero-point. The results 
are reported in Table [T] 

To further check our photometric calibration, in particular for 
stars in fields observed in clear nights, we derive the median stellar 
loci of SDSS and ShaSS stars in (g-r, r—i) colour-colour space (see 
Fig.E). The distance between the two loci, over the whole range of 
colours, never exceeds 0.015 mag, which is less than the expected 
error on colours. 

When photometric calibration is performed, t he bac kground 
is removed with the software SWARP (Bertin et al. |2 oo 3R which 
also standardizes the zero-point to a value of 30.0 mag. 


3 CATALOGUE CONSTRUCTION AND STAR/GALAXY 
CLASSIFICATION 

The procedure adopted for the source extraction is optimized for 
the goals of the survey. Due to the image depth, sources span a 
wide range of size, luminosity and morphology, and thus we need 
a multifaceted approach to obtain robust measures of their aperture 
and total magnitudes. Moreover, the SSC is located at relatively 
low galactic latitude {b ~30), which implies the presence of a large 
number of stars across the survey area, making the star/galaxy clas¬ 
sification a crucial issue for the catalogue’s construction. 

The p hotometric catalogues are produced using the software 
SExtractor jBertin & Amoutslll996ll in conjunction with PSFE}0 
jBertinl I 2 OI 1 1. which performs PSF fitting photometry. We ex¬ 
tract independent catalogues in each band, which are th en matched 
across the four wavebands using STILTS (Tavlor l200^ . 


^ Available at http://www. astromatic .net/so£tware/scamp. 
^ Available at http://www. astromatic .net/so£tware/swarp. 
Available at http://www. astromatic .net/so£tware/ps£ex. 



Figure 3. Colour-colour (g-r,r-i) diagram showing a sample of stars from 
ShaSS (black dots) overimposed on a sample of SDSS stars (red dots). The 
curves trace the loci of the stars in the plane (ShaSS: cyan; SDSS: yellow). 
The distance between the curves never exceeds ~0.015 mag. 


3.1 Catalogue construction 

The procedure adopted for source extraction aims: i) to detect as 
many sources as possible, while minimizing the contribution from 
spurious objects, ii) to produce accurate measurements of positions 
and photometric quantities. Hi) to flag objects in the haloes of bright 
stars and hence could have had their photometry affected. During 
the construction of the catalogues, the results have been always vi¬ 
sually inspected on the images to check the residual presence of 
spurious objects or misclassified objects, like traces of satellites, 
fake objects due to cross-talk, effects of bad columns. 


3.1.1 Source detection 

Sources included in the final catalogue are extracted in four steps: 
i) sky background modelling and subtracting, ii) image filtering. Hi) 
thresholding and image segmentation, iv) merging and/or splitting 
of detections. 

In order to optimize the automatic background estimation, we 
obtained catalogues by adopting different BACK.SIZE and BACK 
FILTERSIZE and compared sources extracted in each catalogue, 
both in terms of number of spurious detections and photometric 
quantities, such as aperture and Kron magnitudes and flux radius. 
Finally, given the average size of the objects, in pixels, in our im¬ 
ages, and in order to minimize the number of spurious detections, 
we set the BACK.SIZE and BACK FILTERSIZE to 256 and 4 respec¬ 
tively, for all fields and bands. To get accurate background values 
for the photometry, the background is also recomputed in an area 
centred around the object in question, setting BACKPHOTO.TYPE 
to LOCAL and the thickness of the background LOCAL annulus 
(BACKPHOTO.THICK) to 24. 

Once the sky background is subtracted, the image must be fil¬ 
tered to detect sources. To determine whether the Gaussian or top- 
hat filter was optimal for our scientific objectives, a specific analy¬ 
sis was carried out on the g-, r- and t-band images for field 8. The 
catalogues produced using the top-hat filter were found to contain 
more spurious sources than those with the Gaussian filter, while the 
photometric measurements as well as the completeness of the cata- 
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Parameter 

Values 

DETECT JIINAREA 

5 

DETECT.THRESH 

1.5(T 

ANALYSIS.THRESH 

1.5(T 

DEBLEND.NTHRESH 

32 

DEBLEND JtINC0NT;,o, 

0.001 

DEBLEND JtINCONTcoW 

0.01 

BACK.SIZE 

256 

BACK.EILTERSIZE 

4 


Parameter 

Values 

PHOT APERTURES 

1.5,3.0,4.0,8.0,16.0,30.0.45.0,90.0 

PHOT APERTURES 

3*FWHM,8*FWHM 

PHOTAUTOPARAMS 

2.5, 3.5 

PHOT_PETROPARAMS 

2.0, 3.5 


Table 3. Adopted aperture diameters, MAGJVUTO and 11AG_PETR0 parame¬ 
ters. PHOT_APERTURES are expressed in arcsec (first row) and as a function 
of the measured FWHM (second row). 


Table 2. Main input parameters set in the SExtractor configuration files. 
DETECT JIINAREA and BACK.SIZE are expressed in pixels. 

logues were confirmed to be equivalent. Finally we chose to apply 
the Gaussian filter (gauss_3. S_5x5 . conv). 

The detection process is mostly controlled hy the threshold¬ 
ing parameters (DETECT.THRESHOLD and ANALYSIS.THRESHOLD). 
The choice of the threshold must be carefully considered. A too 
high threshold results in the loss of a high number of sources in the 
extracted catalogue, while a too low value leads to the detection of 
spurious objects. Hence, a compromise is needed hy setting these 
parameters according to the image characteristics, the background 
rms, and also to the final scientific goal of the analysis. 

For our analysis, the threshold value was chosen to maximize 
the number of detected sources, while simultaneously keeping the 
number of spurious detections to a minimum. To identify the op¬ 
timal threshold value, we counted the number of extracted sources 
for different threshold values in 1 deg^ r-band image and the cor¬ 
responding negative image, which is the scientific image multi¬ 
plied by -1. As the threshold decreases, the number of detected 
sources (real plus spuriuos detections) increases, as shown in the 
left-hand panel of Fig. |4] In the negative image (right panel of 
Fig.® the number of sources (spurious detections in this case) in¬ 
creases smoothly down to a certain value of the threshold, beyond 
which it shows a dramatic change in steepness. The suitable thresh¬ 
old value corresponds to this change in the trend of source number 
counts. As pointed out above, the catalogue was visually inspected 
to avoid residual spurious detections and to verify the deblending 
parameters. 

Two or more very nearby sources can be mistakenly detected 
as a unique connected region of pixels above threshold and, in or¬ 
der to correct for this effect, SExtractor adopts a deblending method 
based on a multi-thresholding process. Each extracted set of con¬ 
nected pixels is re-thresholded at N levels, linearly or exponentially 
spaced between the initial extraction threshold and the peak value. 
Also here we should find a compromise for the deblending param¬ 
eter, since too high a value leads to a lack of separation between 
close sources, while too low a value leads to split extended spiral 
galaxies into several/multiple components. 

In general, the choice of the deblending (but also the thresh¬ 
old) parameter is related to the main scientific interest. If we are 
mostly interested in the analysis of faint and small objects, we have 
to fix low values for deblending and threshold parameters, at the 
cost of splitting up the occasional big bright spiral galaxy into many 
pieces. On the other hand, by fixing larger values for deblending 
and threshold parameters, we correctly deblend the larger objects, 
but we lose depth. However, the present data cover large areas and 
the science involves the analysis of both large and bright as well as 
faint and small galaxies. 

Eor this reason it is not possible to fix a unique value for 
threshold and deblending parameters, and we used a two-step ap¬ 


pr oach in the catalogue extraction (e.g.: Rix et al. 120041 Caldwell et 
al. l2008h . First, we ran SExtractor in a so called cold mode where 
the brightest and extended sources are properly deblended; then, 
in a second step, we set configuration parameters in order to de¬ 
tect fainter objects and to properly split close sources {hot mode). 
We combined the two catalogues by replacing extended objects, 
properly deblended in cold mode, in the catalogue of sources de¬ 
tected in the hot mode, and by deleting multiple detections of these 
extended sources. In order to combine these two catalogues we 
compare the segmentation maps obtained with cold and hot modes. 
Eirst we flagged the different objects detected in the hot mode, but 
lying inside the Kron area of the same object detected in cold mode, 
which are typically 80-120 per field. These objects are, in general, 
large spiral galaxies, galaxies with superimposed smaller galaxies 
or stars and close objects not properly deblended in cold mode. 
Erom the visual inspection of all the flagged objects we chose cold 
mode detection for spiral galaxies and hot mode for other cases. 
The main values set for SExtractor are listed in Tablej^ 


3.1.2 Aperture photometry and total magnitudes 

Measurements of position, geometry, and photometric quantities 
are included in the final catalogues for all the detected and properly 
deblended sources. Among the photometric quantities calculated 
are: aperture magnitudes (MAGJVPER), measured in 10 circular aper¬ 
tures whose fixed diameters are reported in Table|3] isophotal mag¬ 
nitudes (MAG.ISO), computed by considering the threshold value 
as the lowest isophote and using two different elliptical apertures: 
the Kron (HAG AUTO, Kron Il980l) and the Petrosian (HAG_PETR0, 
Petrosian I 19761) magnitudes, which are both estimated through an 
adaptive elliptical aperture. The adopted Kron and Petrosian fac¬ 
tors, and the minimum radius fixed for the elliptical apertures are 
reported in Table |3] (PHOTJIUTOPARABS and PH0T.PETR0PARAMS, 
respectively). 

In order to test the effect of seeing on the ap erture magnitudes 
we used the program SYNMAG (Bundy et ah, l20I2h . With this 
software we can input magnitudes, effective radii and Sersic indices 
to obtain the synthetic magnitude of the source inside an aperture of 
fixed radius. We found that for photometry measured within fixed 
apertures of 8 " diameter, the level of flux lost remains below 0.02 
mag for both exponential (nse„,v; = l) and de Vaucouleurs {nsersic=4-) 
profiles, for seeing levels in the range 0.5-1.1 arcsec and effective 
radii in the range 3.0-100.0 arcsec. 

To describe the size of the sources we extract: the half flux 
radius, i.e. the FLUX_RADIUS containing the 50% of the total light 
and that containing the 90% of the light. We also measured the 
PETRO_RADIUS, defined as the point in the radial light profile 
at which the isophote at that radius is 20% of the average sur¬ 
face brightness within that radius, and the KRON_RADIUS, which 
is the characteristic dimension of the ellipse used to calculate the 
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Figure 4. Number of extracted sources for the r-band image {left-hand panel) and negative image (right-hand panel, see text) of field 12 as a function of the 
SExtractor threshold. The dashed vertical line indicates in both panels the adopted detection threshold. 


Pai'ameter 

Values 

BASIS.TYPE 

PIXELAlUTO 

BASISJJUMBER 

20 

PSE ACCURACY 

0.01 

PSF.SIZE 

25,25 


Table 4. Main input parameters set in the PSFEx configuration file. 
PSF.SIZE are expressed in pixels. 


KRON_MAGNITUDE. The Kron and the Petrosian magnitudes are 
measured within elliptical apertures, whose semi-major and semi¬ 
minor axes are equal to A_IMAGE and B.IMAGE multiplied by the 
KRON_RADIUS and PETRO_RADIUS parameters, respectively. 

As position parameters, we measured the barycenter co¬ 
ordinates both in pixels (X.IMAGE, Y_IMAGE) and in degrees 
(ALPHA.J 2008, DELTA.J 2000), computed as the first order mo¬ 
ments of the intensity profile of the image, and the position of the 
brightest pixel (ALPHAPEAK.J2O0O, DELTAPEAK.J20OO). 

3.1.3 Model-derived magnitudes 

By using the PSFEx tool [BertinI ( 1201 ih it is possible to model the 
PSF of the images. We considered only non-saturated point sources 
with a SNR higher than 80. Spatial variations of the PSF were mod¬ 
elled with a third degree polynomial of the pixel coordinates. The 
main values set for PSFEx parameters are given in Table |4l More 
details on th e PSF model l ing w ith the PSFEx t ool can be found 
[Bertiij j201lh . lMohr et alJ ( 1201 2h and Bouy et al. ( |2013|) . 

Hence, it is possible to run SExtractor taking the PSF models 
as input, using them to carry out the PSF-corrected model fitting 
photometry for all sources in the image. With a combination 
of SExtractor and PSEEx, we obtained magnitudes from: (i) 
the PSF fitting (MAG.PSF and the point source total magnitude 
MAG.POINTSOURCE); (ii) the fit of a spheroidal component 
(MAG.SPHEROID); (Hi) the fit of a disc component (MAGJIISK); and 
(iv) the sum of the bulge and the disc components, centred on the 
same position, convolved with the local PSF model (MAG.M0DEL). 
We also extracted morphological parameters of the galaxies, 
such as spheroid effective radius (SPHEROID.REFE.IMAGE), 
spheroid Sersic index (SPHEROID.SERSICN), spheroid 


aspect (SPHEROID.ASPECT.IMAGE), disc scale length 
(DISK.SCALE.IMAGE) and disc aspect (DISKJISPECT.IMAGE). 

The model of the PSF is also helpful to get a more accurate 
star/galaxy classification (see Sect. l3.3b . 


3.2 Haloes and spikes of bright stars, and image borders 

Multiple reflections in the internal optics of OmegaCAM can pro¬ 
duce complex image rings and ghosts (hereafter star haloes) near 
bright stars These haloes are characterized by a central region 
in which the surface brightness is depressed and an outer corona 
with enhanced surface brightness, both characterized by an irreg¬ 
ular pattern as shown in Fig. [5] (left panel) for the i-band images. 
Those haloes all have the same radius of 830 pixels in all bands 
and the position of their centres is related to the position of the 
parent star in the OmegaCAM field. The ‘intensity’ (depth of the 
inner depression and brightness of the corona) of the halo depends 
on the band, being more manifest in the i band and practically neg¬ 
ligible in u band. Moreover, the intensity of the halo is related to 
the magnitude of the parent star. To quantify this behaviour we as¬ 
signed a value proportional to the halo masks: mask=INT(100x 
imag)- Notice that parent stars are saturated in our images, so we 
draw their magnitudes from the USNO-Bl.O catalogue. Whenever 
two or more haloes overlap, the mask assumes locally the value 
corresponding to the halo of the brightest star. We limited the cre¬ 
ation of masks to stars brighter than i=9.6 mag, although the haloes 
produced by stars with i > 9.0 mag already do not show surface 
brightness excesses higher than ~0.2 times the r.m.s. of the back¬ 
ground. In all 23 VST fields in i band, we identified and masked 
299 haloes. In the right panel of Fig. [5] we show an example of a 
halo mask. 

The saturated stars with haloes also affect the images by pro¬ 
ducing spike features which are masked as follows. First we iden¬ 
tified the position of all saturated pixels according to the value cor¬ 
responding to the SATURATE keyword in the header of each image. 
Then we applied a region growth algorithm around each saturated 
pixel to obtain a mask consisting of three levels, corresponding to 


^ See http://www.eso.org/observing/d£o/quality/OMEGACAM/qc/problems.html 
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Figure 5. A crowding of five saturated stars with related haloes (left) and 
the corresponding halo mask (right). The i-band magnitudes of the stars are, 
left to right, 8.77, 8.24, 8.99, 8.56, 7.81. A logarithmic stretch was applied 
to the image to increase the background. The image corresponds to an area 
of 9' X 6'. 



* 

IP 


Figure 6. An area containing a star with prominent spikes and other sat¬ 
urated stars (top) and the comesponding spike mask (bottom). The three 
levels of the mask correspond to three different surface brightness thresh¬ 
olds (/i/ < 21, 20.5, 20 mag arcsec^^, coded here in blue, green and red 
respectively). The shown field measures 2.8' X 7.7'. 

three threshold values of the surface brightness of 20, 20.5 and 21 
mag arcsec^^. An example of a spike mask is shown in Fig.|^ 

The halo and spike masks have been used to flag each source 
s, considering the circular area embedding 50% of the flux, A50, . 
From both halo and spike masks we derive a couple of flags: i) the 
fraction of A50, affected by the halo and/or spike mask; ii) the halo 
and/or mask value within the affected portion of A50,. 

The dither acquisition mode adopted for the observations pro¬ 
duces a non-uniform coverage of the co-added image, so that the 
catalogue of detections covering the whole mosaic image will con¬ 
tain detections from regions with different depths. For a five(nine)- 
dither pattern at least three(five) exposures cover almost the whole 
image in ugi(r). Along the borders, however, there are stripes where 
only 1-2(1-4) exposures contribute to the final co-added mosaic in 
ugi(r). We decided to exclude these detections from the final cata¬ 
logues. This cut was possible thanks to the 3 overlaps among the 
contiguous 1 deg^ VST fields, which ensure that a common area 
(and common detections) is present even after the borders are cut. 
To mask out we empirically identified the external areas with dif¬ 
ferent exposures by just associating a range of values of the weight 



Figure 7. Cutoff of an exposure mask near to a comer of a VST image. 
The area shown is about 14' X 12'. The colours correspond to the different 
numbers of exposure as indicated. 

mask to a number of exposures. An example of an exposure mask is 
given in Fig.|7] Objects lying in areas with mask values lower than 
3(5) are not included in the survey catalogues, except for a very few 
number of sources located in the gaps among the detectors. 

We also excluded those objects for which more than 80% of 
their pixels (within a circular region with radius equal to 50% flux 
radius) have zero weight. 

3.3 Star/Galaxy separation 

To separate extended and point-like sources, we adopt a progres¬ 
sive approach analogous to those described in Annunziatella et 
al. ( I 2 OI 3 I) , using the following parameters provided by SExtrac- 
tor: (i) the stellarity index (CLASS_STAR); (ii) the half-light radius 
(FLUX_RADIUS); (iii) the new SExtractor classifier SPREAD_MODEL; 
(iv) the peak of the surface brightness above background f/imax); (v) 
a final visual inspection for objects classified as galaxies but with 
borderline values of the stellarity index (CLASS_STAR>0.9; see be¬ 
low). 

The stellarity index results from a supervised neural network 
that is trained to perform a star/galaxy classification. It assumes 
values between 0 and 1. In theory, SExtractor considers objects 
with CLASS.STAR equal to zero to be galaxies, and those with value 
equal to one as stars. As an example in Eig.[8^ CLASS_STAR is plot¬ 
ted as a function of the Kron magnitude for r-band sources in field 
12. The sequence of unsaturated stars (r > 15.5 mag) is clearly sep¬ 
arated from galaxies by selecting a CLASS.STAR value above 0.98 
only down to r=20.5 mag. For magnitudes fainter than this value, 
lowering the established limit to separate stars and galaxies causes 
an increase of the contamination of the star subsample from galax¬ 
ies. So, we can only adopt this parameter to classify bright sources 
(r <20.5 mag) as stars, when CLASS_STAR>0.98. 

By using (FLUX_RADIUS) as a measure of source concentra¬ 
tion, we can extend the classification to fainter magnitudes. Fig¬ 
ure]^ shows that the locus of stars, defined according to the rela¬ 
tion between half-light radius and Kron magnitude is recognizable 
down to /■=22.8 mag. This limit is 0.7 mag brighter than the com¬ 
pleteness limit of the r-band catalogue (see Sect. l4.3t . For this rea¬ 
son we used the new SExtractor classifier, SPREAD_MODEL, which 
takes into account the difference between the mod el of the source 
and the model of the local PSF foesai et ani2012l) . to obtain a re- 






















8 A. Mercurio et al. 


liable star/galaxy classification for the faintest objects in our cat¬ 
alogue. By construction, SPREAD_MODEL is close to zero for point 
sources, positive for extended sources (galaxies), and negative for 
detections smaller than the PSF, such as cosmic rays. Figure]^ 
shows the distrihution of the SPREAD_MODEL as a function of Kron 
magnitude. Stars and galaxies tend to arrange themselves in two 
different places of the plot distinguishable up to t~23.8 that is a 
limit 0.3 mag fainter than the completeness limit of the catalogue. 
Based on this diagram, we classified as galaxies all sources with 
SPREADJVIODEL > 0.005 and r <23.8, or SPREADJVIODEL > 
0.003 and r >23.8. 

Einally, in Eig.[^we plot yu^ax as a function of the Kron mag¬ 
nitude. This plot is used in order to select saturated stars (vertical 
dashed line in all panels of Eig. 8). 

Since the r hand is the deepest band of the survey, and the one 
conducted in the best seeing conditions, it will be used for clas¬ 
sification of sources in the cross-correlated catalogue of the four 
bands. The same criteria have also been used for the i and g band 
independently and the results are consistent. 

The star/galaxy separation of the u hand relies on the r hand, 
since the nature of the w-hand emission, which is very sensitive to 
the presence of SE regions. Eor this band a careful cross-correlation 
with the other catalogues is also required since the dehlending pa¬ 
rameter cannot be suitably calibrated. We will further detail the 
construction of the K-band catalogue in a dedicated article. 

After completing the star/galaxy separation, a visual in¬ 
spection of those objects classified as galaxies but with 
CLASS_STAR>0.9 (which are about 10 per field), is performed. 

To check the reliability of the star/galaxy classification in the 
i-band images, we performed simulations by adding artificial stars 
and galaxies, through a stepwise procedure. To preserve the over¬ 
all source density, we split the artificial stars and galaxies into 
eight magnitude bins of width 0.5 mag, over the magnitude range 
18< i <22, producing eight different simulated images for each real 
image. To define the PSE we take advantage of the PSE modelled 
according to PSEEx. To simulate galaxies with dimensions repre¬ 
sentative of real deep images, we defined the size in pixels of the 
model in the output image, according to the distribution of the mea¬ 
sured half-flux radius as a function of the total magnitude obtained 
for the real galaxies. We attempted to recover artificial sources by 
running SExtractor again with the same parameters for object de¬ 
tection and classification as on the original images. We were able 
to correctly separate 100% of simulated stars and galaxies down 
to 1=20.0 mag. We erroneously classify as galaxies 0.8%, 1.5%, 
2.1% and 3.0% of stars in the magnitude bins 20.0-20.5, 20.5-21.0, 
21.0-21.5 and 21.5-22.0, respectively. All galaxies were correctly 
recovered down to 1=21.0, while in the last two magnitude bins just 
0.7% and 2.2% were misclassified as stars. 

To test if the method adopted to classify stars and galaxies 
could introduce a bias in the distribution of galaxy sizes we com¬ 
pared our observed distribution of real galaxy sizes, with those de¬ 
rived by using completely independent methods. In particular, we 
verified that our size distribution for galax ies dow n to M*-l-6 was in 
agreement with those of D’Onofrio et al. ( l2014h for WINGS clus¬ 
ters. 


However, compact galaxies might be confused with stars due 
to their small radii and comparably high surface brightness. Pre¬ 
vious studies on the sizes of dwa rf and compact d warf g alaxies 
in C oma (e .g. Graham & Guzman l2003l . Price et al. 120091 Hoyos 
et al. [20^ and in other nearby clusters li ke Cen taurus (Misgeld 
et al. 20091) . Antlia (Smith Castelli et al. |2012|) . Virgo, Eomax 
and Perseus (Weinmann et al. 1201 ih have shown that the effective 


Table 5. ShaSS photometry. 


Band 

Exp. Time 
[s] 

completeness 

[mag] 

detection [5 (t] 
[mag] 

seeing 

[arcsec] 

u 

2955 

23.8 

24.4 

0.9 

g 

1400 

23.8 

24.6 

0.7 

r 

2664 

23.5 

24.1 

0.6 

i 

1000 

22.0 

23.3 

0.7 


a) Median FWHM estimated from the already obseived fields. 


radii of these galaxies can vary from few kiloparsec to hundreds of 
parsec. To investigate to what extent we are able to properly dis¬ 
tinguish compact galaxies from stars we run ad hoc simulations. 
We added to the images artificial compact galaxies and used the 
same detection method as for the real catalogues. Simulated galax¬ 
ies with a minimum effective radius of 0.7 kpc (i.e. 1 arcsec at 

z=0.048) are all recovered. Then, considering galaxies with a min¬ 
imum effective radius of 0.6 kpc, we erroneously classify as stars 
only 2.2% of compact galaxies. Einally, by simulating galaxies with 
effective radii down to 0.4 kpc we were able to classify as extended 
objects 60-70% of galaxies. However, we would like to underline 
that 400 pc corresponds to ~ 0.4 arcsec at the mean distance of the 
SSC (z=0.048), so, to assess the nature (star or galaxy) of these ob¬ 
jects, we need either a measure of the redshift or to perform a fit to 
the surface brightness distribution. 


4 PHOTOMETRY 

In this section, we will focus on the photometric depth, the accu¬ 
racy of the derived photometry and the completeness magnitude. 
All the quantities shown in this section are obtained by excluding 
saturated stars, sources alfected by saturation spikes, stars in haloes 
of bright stars, and those with less than half the total exposure of 
the frames, which means that we included the sources with at least 
three exposures in the ugi bands and with at least five exposures in 
the r band. We remind the reader that the results for the u band are 
preliminary. 

The results of Sections [4.1 1 and [43] are based on the ugri cat¬ 
alogues including 11, 10, 14 and 23 fields, respectively. Eor each 
band we cross-correlate the field catalogues with STILTS selecting 
in the overlap regions the detection from the image with the best 
seeing. This criterion is also adopted for producing the final i-band 
catalogue covering the whole ShaSS area. 

We also report in Table |5] a summary of the exposure time, 
median seeing, completeness limits and detection limits at 5cr for 
each band. 


4.1 Photometric depth 

We estimated the photometric depth of each field by randomly plac¬ 
ing 10000 3-" apertures on the VST image of each field. Prom 
the resulting standard deviation in the flux measurements obtained 
within these apertures, the corresponding 2Qa and 5cr magnitude 
limits were determined for each image. The mean values obtained 
in each field are 22.9, 23.1, 22.6, 21.7 in ugri band at 20cr and 
24.4, 24.6, 24.1, 23.3 in ugri band at 5cr. The variation of these de¬ 
tection limits across the fields is less than 0. 05 in ug, and 0.1, 0.2 in 
r and i band, respectively. Merluzzi et al. ( 1201 5h reported slightly 
different values since they measured the SNR within 3 " diameter 
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Figure 8. Two-dimensional histogram of SExtractor stellarity index (panel a), half-light radius (panel b), spread model (panel c), and Pmax (panel d) as a 
function of Kron magnitude for sources in r-band image of field 12. The vertical dashed line in all panels is the magnitude limit adopted for saturated objects. 
The horizontal lines in panel c) mark the limit used to separate stars and galaxies (see text). The points are colour coded according to their number counts as 
reported in the vertical colour bars. 


apertures as a function of the magnitude, according to equation (61) 
of SExtractor User’s Manual, while in this paper we measured the 
detection limits directly on the images. 


4.2 Photometric accuracy 

To test the accuracy of our photometry we used the stars in the 
3 arcmin wide stripes of overlap between adjacent VST fields. The 
magnitudes adopted for this comparison are those derived in an 8" 
diameter aperture, which ensures that our estimates are not affected 
by the seeing differences between fields. The average number of 
stars in each strip are ~650 in u and 1500 in gri bands. 

We computed the magnitude difference Aj,^ =m^- m^ as a func¬ 
tion of magj = |x(m^-l-m^) for all stars belonging to each pair ij of 
adjacent frames. Here i and j identify the different VST fields and 
k refers to the stars, so that m^ is the magnitude of the A;-th star in 


field i and is the difference in magnitude of the A:-th star between 
fields i and i. 

To estimate the uncertainties on magnitudes, we collected the 
absolute values of the differences |A^''| and mag^ of all stripes into 
two single vectors A^ and mag*. We then divided the sample into 
N (-26) equally populated ^ 1000-2000) bins of magnitude, 

and for each bin, we computed the standard deviation of A*,, ctj,,,, 
using a 3cr rejection for the outliers (the subscript bin refers to the 
different bins in which the sample was divided). 

We adopt cr^,,, as the measure of the uncertainties of the mag¬ 
nitude differences in each bin, and therefore 1 / V2 x cri,„ as an 
empirical measure of the uncertainties of the magnitudes. These 
uncertainties are shown in the upper panel of Fig. as functions 
of magj. The ratios of the empirical uncertainties with those es¬ 
timated by SExtractor are shown in the lower panel of the same 
figure. 

The ratios of the two uncertainties may be expressed as a func¬ 
tion of magnitude and band as 

^111 1 Extractor — 'X eXp(bj^ X ITl) + Cj^ , 


® https: //WWW. astromatic.net/pubsvn/so£tware/sextractor/trunk/do(l>fe'extrfalztt,cnj,p((M= 1-4) are coefficients depending on the band 
















10 A. Mercurio et al. 



j_^ ^ _I_^_I_ ^ _I_ I _ ^ _I_I_ ^^^ _ s_a 
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Magnitude 


Figure 9. Top\ magnitude uncertainties derived from stars in the 3 arcmin 
overlapping stripes between adjacent VST fields (continuous lines) and un¬ 
certainties estimated by SExtractor (dashed) as functions of the magnitude. 
The colours code the wavebands as indicated in the lower panel. Bottom: 
ratio between the two uncertainty estimates. Continuous lines are the fit to 
the real data (filled circles) whose coefficients are reported in Tab[6] 


Table 6. Coefficients of the ratio between empirical and SExtractor uncer¬ 
tainties as a function of the magnitude - R=aiXexp[a 2 X (mag-20)]-tas. 


Band 

ai 

a2 

as 

u 

1.26 

-0.985 

1.353 

g 

18.55 

-0.762 

1.144 

r 

9.39 

-0.839 

1.397 

i 

3.67 

-0.858 

1.467 


which were derived by least-square fits to the curves in Fig.[9](bot- 
tom). 

The uncertainties given in the catalogue are those given by 
SExtractor (i.e. the nominal error l^sExtmctor)- To obtain a more real¬ 
istic error, the nominal error should be multiplied by A„, I As Extractor, 
with the zero-point uncertainty added then in quadrature. In Tab.|^ 
we give the value of the multiplicative factor as a function of the 
magnitude. Zero-point errors are reported in Table[T] 

4.3 Catalogue completeness 

Following the method of iGarilli, Maccagni & AndreonI l ll999l) . we 
estimated the completeness magnitude limit as the magnitude at 
which we begin to lose galaxies because they are fainter than the 
brightness threshold inside a detection aperture of 1.5 " diametej^, 
for the g, r and i bands. For the u band we give a first estimate 
of the completeness, based on the distribution of the Kron magni¬ 
tude, since the analysis of the photometry is still ongoing. In all 
panels of Fig. [la the vertical blue dashed lines represent the de¬ 
tection limit, while the red continuous lines are the linear empiri- 

^ This aperture was adopted being suitable for all images and comparable 
to the larger seeing value. 



Magnitude inside 1.5 arcsec diameter 

Figure 10. Distribution of the SExtractor magnitude inside a 8.0" diameter 
as a function of the magnitude inside a detection aperture of 1.5" diameter 
for gri bands. The horizontal and vertical blue dashed lines indicate the de¬ 
tection and the completeness limits, respectively. The red continuous lines 
are the linear relation between the magnitude within 8.0" aperture diameter 
and the magnitude within the detection aperture, minus/plus Icr (red dashed 
lines). 


cal relation between the magnitude within 8.0" diameter aperture 
and the magnitude within the detection aperture. The relation be¬ 
tween the two magnitudes shows a scatter, depending essentially 
on the galaxy profiles. Taking into account this scatter (see dashed 
red lines in Fig. llOt . we fixed as a completeness magnitude limit 
(blue dashed horizontal line) the intersection between the lower Itr 
limit of the relation and the detection limit, which corresponds to 
23.8, 23.8, 23.5, 22.0 in ugri bands, respectively. We checked the 
percentage of galaxies retrieved at the completeness limit defined 
above, by means of the recovery rate of artificial galaxies inserted 
in the images and retrieved with identical procedure as those used 
for real sources. We simulated 100 galaxies for each 0.5 mag bin 
over the magnitude range 18< i <23.5 mag and effective radii rang¬ 
ing from 0.3 to 20 arcsec. We verify that the catalogues turned out 
to be 93% complete at the total magnitudes of 23.8, 23.8, 23.5, 22.0 
in ugri band, respectively, and they reach the 89% of completeness 
0.5 mag fainter. 

Fig.E] shows the distribution of the Kron magnitude, where 
vertical continuous lines indicate the limit to which the catalogues 
are complete. 
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Figure 11. Number counts of galaxies over the observed fields (11, 10, 14, 
23, for M, g, r and i, respectively) per 0.1 mag bin. Vertical lines mark the 
completeness magnitudes, u- and g-band catalogues have the same com¬ 
pleteness magnitude indicated by the blue vertical line. 


5 DESCRIPTION OF THE DATA BASE 

The (-band catalogue is published and it will be regularly updated, 
through the use of the Virtual Observatory (VO) tools. 

The data model of the ShaSS data base is based on the 
standard Entity-Relationship (ER) paradigm. Its architectural di¬ 
agram is shown in Fig. [m Its physical model (band-specific and 
correlation tables) is instantiated in the logical DataBase Man¬ 
agement System (DBMS), based on the open source mySQL 
data base and on the SVOCat service, a VO software tool de¬ 
veloped in the framework of the Spanish Virtual Observatory 
(http://svo2.cab.inta-csic.es/vocats/SVOCat-doc/) . 

The DBMS provides direct access to data through SQL- 
based queries and VO standard interoperability (i.e. direct ta¬ 
ble/image exchanging among VO tools, such as Topcat and Al- 
adin). The service is publicly accessible via browser at the address 
http://shass.na.astro.it . 

The Shass data base is also publicly available within 
the EURO-VO registry framework, under the INAF- 
DAME Astronomical Archive identification authority 
(ivo: //dame. astro. it/shass- i i. 

Each object in the data base has a unique primary key 
(SHASS.ID), which univocally identifies sources. It is composed 
by a string containing 14 characters, where five characters are for 
ShaSS and none digits NNNNNNNNN represent the internal object 
identification number. 

The SHASS.ID field represents the primary identifier of an ob¬ 
ject within the corresponding band table. Whenever an object has a 
counterpart in different bands, the sequence of different SHASS_IDs 
for each reference band can be obtained by simple queries to the 
correlation table in the data base (for instance the corrjab table in 
Fig.[l2}. In this table the fields group and group.size are reported 
for each band. These indicate the possible presence of a group of 
objects matching with the single source within 1" distance. To each 
group of sources matching with one specific object is assigned a 


unique integer, recorded in the group field, and the size of each 
group is recorded in the group_size field. Sources which don’t 
match any others (singles) have null values in both these fields. 

In the data base we report the barycenter coordinates in de¬ 
grees (RAdeg, DECdeg), the geometrical parameters of the ellipse 
that describes the shape of the objects: semi-major and semi-minor 
axes, and position angle (A, B, THETA). As an indicator of the area 
covered by the objects we include the half flux radius (FR 50 ) which 
corresponds to the radius of the isophote containing half of the 
total flux. We add the Kron radius (RK) and the Petrosian radius 
(RPETRO), which are the indicators of size of the Kron and the Pet¬ 
rosian aperture respectively. Both of these radii are expressed in 
multiples of the semi-major axis. 

Among the extracted photometric quantities, we report in the 
public catalogue three aperture magnitudes (MA) (see Table |7j, the 
Kron magnitude (MK), the magnitude resulting from the PSF fitting 
(MPSF), the model magnitude obtained from the sum of the spheroid 
and disk components of the fitting (MMODEL), and the Petrosian 
magnitude (MPETRO). Magnitudes are not corrected for galactic ex¬ 
tinction and we give the relative uncertainties as derived by SEx- 
tractor (but see also Sect. l4.2t . 

Finally we provide the stellarity index from SExtractor (SI) 
and five flags: the star/galaxy flag (SG), the halo fraction flag (HFF), 
the halo flag (HF), the spike fraction flag (SFF) and the spike flag 
(SF). 

The star/galaxy flag is fixed according to the star/galaxy sepa¬ 
ration described in Sect. 13.31 Stars classified according to the SEx¬ 
tractor stellarity index, half flux radius and the spread model have 
SG= 1,2,3 respectively, while those classified through the visual in¬ 
spection have SG=7. Saturated stars are indicated by SG=9. Some¬ 
times values of the star/galaxy flag SG=4,5 could be present. They 
indicate stars aligned in a ‘secondary sequence’ visible in a few 
observed fields in the plots of the half light radius and the ji^ax ver¬ 
sus the Kron magnitude, respectively. The origin of this secondary 
sequence is unclear, but it seems a ’’random effect” since there is 
no correlation with observing night, airmass or spatial position of 
stars, and only constitutes a small number of stars. 

As stated in Sect. l3.3l the r band will be used for the classifica¬ 
tion of sources in the cross-correlated catalogue, while in the g and 
i single-band catalogues we report the classification independently 
derived for each band, which is mostly consistent with that of the r 
band. 

The halo and the spike fraction flag indicate the fraction of the 
object area, defined as the circular area with radius equal to half flux 
radius, affected by the presence of the halo or the spike of a bright 
saturated star. Since, as pointed out in Sect. l3.2l the ‘intensity’ of the 
halo is related to the magnitude of the parent star, flag HF is equal 
to lOOx i„„^, where i„„j is the star magnitude from the USNO-B 1.0 
catalogue. Spikes flag values are equal to 1,2 or 3, corresponding 
to three threshold values of the surface brightness of 20, 20.5 and 
21 mag arcsec^^ of saturated regions (see Sect. l3.2L 

A description of the units used for each of the quantities in the 
data base is reported in Table |7] 

We plan to complement the present data base with the cata¬ 
logues of the M, g and r bands and the cross-correlated catalogue 
following the completion of the observations. 

6 SUMMARY 

The ShaSS will map an area of ~23 deg^ (■'- 260 h/^ Mpc^ at 
z=0.048) of the SSC with the principal aim being to quantify the 
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Parameter 

Units 

Description 

ID 


ShaSS identification 

RAdeg 

deg 

Right ascension (J2000) 

DECdeg 

deg 

Declination (J2000) 

MK 

mag 

Kron magnitude 

EMK 

mag 

Error on Kron magnitude 

RK 


Kron radius expressed in multiples of major axis 

A 

deg 

Major axis 
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deg 

Minor axis 
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Position angle (CCW/x) 
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Aperture magnitude inside 1.5 " diameter 
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Aperture magnitude inside 4.0 " diameter 
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Error on aperture magnitude inside 4.0 " diameter 
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mag 

Aperture magnitude inside 8.0" diameter 
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mag 

En'or on aperture magnitude inside 8.0 " diameter 
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pixel 

Radius of the isophote containing half of the total 
flux 

SI 


SExtractor stellaiity index 

MPSF 

mag 

Magnitude resulting from the PSF fitting 
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En'or on magnitude resulting from the PSF fitting 

MMODEL 

mag 

Magnitude resulting from the model of the spheroid 
and disc components 

EMMODEL 
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En'or on magnitude resulting from the model of the 
spheroid and disc components 

MPETRO 

mag 

Petrosian magnitude 
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mag 

Error on Petrosian magnitude 

RPETRO 


Petrosian radius expressed in multiples of major 
axis 

SG 


Star Galaxy separation 

HEE 


Halo fraction flag 

HE 


Halo flag value 
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Spike fraction flag 

SF 


Spike flag value 


Table 7. Parameters reported in the data base. 


influence of hierarchical mass assembly on galaxy evolution, and 
to follow this evolution from filaments to cluster cores. 

ShaSS will provide the first homogeneous multi-band imag¬ 
ing covering the central region of the supercluster, including optical 
(ugri) and NIR {K) imaging acquired with VST and VISTA, which 
allows accurate multiband photometry to be obtained for the galaxy 
population down to m*-l-6 at the supercluster redshift. In particu¬ 
lar, the r-band images are collected with a median seeing equal to 
0.6 arcsec, corresponding to 0.56 kpc h^Q at z=0.048 and thus en¬ 
abling the internal properties of supercluster galaxies to be studied 
and distinguish the impact of the environment on their evolution. 

In this article, we described the methodology for producing 
photometric catalogues. The optical survey is ongoing and the anal¬ 
ysis presented in this work is performed on 11, 10, 14 and VST 
fields out of 23 in ugr bands, nevertheless in all bands the consid¬ 
ered subsamples are representative of the final whole sample. The 
i-band imaging is instead complete. 

The catalogues are produced using the softw are SExtracto r 
dSertin & Arnoutsiri996h in conjuction with PSFEx (lBertinll201lh . 
and a careful analysis of the software outputs. 

We were able to obtain a robust separation between stars and 
galaxies up to the completeness limit of the optical data, through a 
progressive approach using: i) the stellarity index (CLASS_STAR); ii) 
the half-light radius (FLUX_RADIUS) ; iii) the new SExtractor clas¬ 
sifier SPREAD_M0DEL; iv) the peak of the surface brightness above 


background (/^max); v) a final visual inspection for objects classified 
as galaxies but with CLASS_STAR>0.90. 

The ShaSS catalogues reach average 5cr limiting magnitudes 
inside a 3" aperture of Mgn'=[24.4,24.6,24.1,23.3] and a complete¬ 
ness limit of Mgn=[23.8,23.8,23.5,22.0], which corresponds to 
m*r+8.5 at the supercluster redshift. These values correspond to 
the survey expectations. 

The i-band catalogue is released to the community through the 
use of the Virtual Observatory tools. It includes 734 319 sources 
down to i-22.0 over the whole area. The catalogue is 93% com¬ 
plete at this magnitude limit and 34% of the sources are galax¬ 
ies. The service is publicly accessible via browser at the address 
http://shass.na.astro.it . The Shass data base is also pub¬ 
licly available within the EURO-VO registry framework, under 
the INAF-DAME Astronomical Archive identification authority 
(ivo://dame.astro.it/shass-ii. 
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